Memory effects in attenuation and amplification quantum processes 
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With increasing communication rates via quantum channels, memory effects become unavoidable whenever 
the use rate of the channel is comparable with the typical relaxation time of the channel environment. We then 
introduce a model of bosonic memory channel, describing correlated noise effects in quantum optical processes 
via attenuating or amplifying media. To study such a channel model we make use of a proper set of collective 
field variables, which allows us to unravel the memory effects, mapping the n-fold concatenation of the memory 
channel to a, unitarily equivalent, direct product of n single-mode bosonic channels. We hence estimate the 
channel capacities by relying on known results for the memoryless setting. Our findings show that the model 
is characterized by two different regimes, in which the cross-correlations induced by the noise among different 
channel uses are either exponentially enhanced or exponentially reduced. 

PACS numbers: 03.67.Hk, ()5.4().Ca, 42.50.-p, 89.70.-a 



I. INTRODUCTION 

Any quantum process can be described as a completely pos- 
itive and trace preserving map (quantum channel) acting on 
the set of trace class operators (quantum states). In quan- 
tum information theory, one of the most important problems is 
finding the maximal rates (i.e., the capacities) at which quan- 
tum or classical information can be transmitted via quantum 
channels with vanishingerror in the limit of large number of 
signals (channel uses) 111]. Earlier works on the subject fo- 
cused on models where the noise affecting the communication 
is assumed to act independently and identically for each chan- 
nel use (memoryless quantum channels). Recently, however, 
an increasing attention has been devoted to correlated noise 
models (quantum memory channels), see, e.g., jj] and ref- 
erences therein. Memory effects in the communication may 
arise when the dynamics of the media is characterized by 
temporal correlations which extend on timescales which are 
comparable with the time delay between consecutive chan- 
nel uses. For instance, under certain conditions optical fibers 
may show relaxation times or birefringence fluctuations times 
longer than the separation between successive light pulses ytl- 
Similar effects occur in solid state implementations of quan- 
tum hardware, where memory effects due to low-frequency 

Further- 
memory 



impurity noise produce substantial dephasing | 
more, moving from the model introduced in 
noise effects have also been studied in the contest of many- 
body quantum systems by relating their properties to the cor- 
relations of the channel environmental state Oa] or by studying 
the information flow in spin networks 113]. 

It is generally believed that memory effects should im- 
prove the information transfer of a communication line ]J8|]. 
However, finding optimal encodings is rather complex and 
up to date only a limited number of models have been 
solved ]|2, S UOtl . In the context of bosonic channels lUlll . 
we have recently introduced an exactly solvable model which 
effectively describes the transmission of signals along attenu- 
ating media characterized by finite relaxation times. It is char- 



acterized by two parameters ^, k G [0, 1], enabling to describe 
different scenarios ranging from memoryless to inter-symbol 
interference channels Ill2ll . up to perfect memory configura- 
tions II13I1 . providing thus the first comprehensive quantum in- 
formation characterization of memory effects in the continu- 
ous variable setup. 

In this article we show the details of such a model and we 
extend it to encompass amplification processes besides attenu- 
ation ones (this is formally obtained by extending the parame- 
ter space to include also k > 1 values). For such processes we 
estimate the classical and the quantum capacity Ill4ll showing 
their enhancement with respect to the memoryless case and 
proving that, for the attenuating memory channel, coherent 
state encoding is optimal for transmitting classical informa- 
tion. This is accomplished by applying suitably chosen en- 
coding and decoding unitary transformations which allow us 
to unravel the correlations among different channel uses. In 
such a way the channel is mapped into a unitary equivalent 
one, in which noises affecting different channel uses are in- 
dependent, although not identical II15I1 . Using this mapping, 
the channel capacities can be estimated by relying on known 
results on memoryless channels, which in the limit of large 
number of channel uses provide converging lower and upper 
bounds. Our results also allows us to point out the existence 
of two different regimes characterized by a threshold value of 
the product /.ik. Indeed, even in the presence of amplification, 
for UK < 1 the correlations among different channels uses in- 
duced by the channel dynamics are depressed exponentially 
with the distance among the uses. This strongly resemble the 
behavior of the so called forgetful memory channels, which 
were introduce in Ref. ^ in the context of finite dimensional 
systems. Vice-versa for /ik > 1 the memory effects, as well 
as the correlations between different channel uses, are expo- 
nentially enhanced by the amplification. Interestingly enough, 
irrespectively of these rather different behaviors the capacities 
of the model show no discontinuities at the critical threshold. 

Even if the model we analyze is presented in an abstract 
setting, we believe that the results presented in the paper, as 



typical in optical communication analysis ||16ll, may have a 
potential impact on a variety of different realistic setups which 
are currently in use, ranging from fiber optics lUTIl . to wave- 
guide |18y, up to free-space signaling I" 19]. Moreover, the 
model can be applied for describing memory effects at the 
interface between matter and light, e.g., in quantum repeaters 
and quantum memories based on ensemble of atoms, in which 
the state of the atomic ensemble is described in terms of an ef- 
fective bosonic degree of freedom |20|1. 

The article is organized as follows. In Sec.|II]we present the 
memory channel model. In Sec.|III]we introduce the technique 
of memory unraveling. In Sec.|IV]we discuss the asymptotic 
properties of the channel in the limit of many channel uses and 
characterize the threshold. The memory unraveling technique 
and the asymptotic properties of the channel are hence used 
for evaluating the classical and quantum capacity in Sec. |V] 
In Sec.|Vl]we argue on the forgetfulness of the model. Finally, 
Sec. lVIII we present concluding remarks. 



II. THE MODEL 

The most relevant effects for optical communications are 
the attenuation and amplification processes, which are mod- 
eled as an exchange interaction of each ingoing mode (sig- 
nal) with a corresponding external bosonic mode (idler), mod- 
eling a local environment. We define a model of Gaus- 
sian memory channel by concatenating these basic transfor- 
mations. An attenuation process involving a pair of ingo- 
ing modes, described through the ladder operators {a,a^} 
and {c, c^}, and the corresponding outgoing modes {b,b^}, 
{d, d^}, is described by a unitary transformation of the form 
Ubs = cxp {9a^c — h.c), modeling a beam-splitter mixing 
the two modes. In the following we assume the parameter 9 
to be real and positive, yielding the Heisenberg-picture trans- 
formations 



b ~ C/bs a Ubs = cos 9a — sin 9 c , 
d = C/bs c Ubs = cos 9c + sin9a, 



(la) 
(lb) 



together with their hermitian conjugates. Different choices for 
the parameter 9 do not affect the main results of our analysis. 
The quantity cos^ 9 is usually referred to as the beam-splitter 
transmissivity. Analogously, we will consider the linear para- 
metric amplifier, described by a unitary operator of the form 
L/la ~ cxp {xo.^c' ~ h.c). For simplicity we assume the pa- 
rameter X to be real and positive, even though other choices do 
not substantially alter our analysis. In the Heisenberg picture 
such a transformation acts on the ladder operators as follows 



b = U-lA a, [/la = cosh x a + sinh x c 
d — UhA c [/la = cosli X c + sinh x a 



(2a) 
(2b) 



together with the hermitian conjugates. The quantity cosh^ x 
denotes the gain of the linear amplifier The circuital repre- 
sentations of these basic unitary transformations are in Fig.[Tl 
To define the action of the quantum memory channel, we 
consider a sequence of n consecutive channel uses. Such a 
sequence is associated with a collection of n bosonic modes 



1^ 




FIG. 1: (Color online.) Basic building blocks of the Gaussian mem- 
ory channel, mixing the signal (represented by the horizontal lines) 
and the idler mode (vertical lines). On the left: the circuital represen- 
tation of the beam-splitter transformations in Eq.s l[T). On the right: 
the circuital representation of the linear amplifier in Eq.s ^. 



m 



FIG. 2: (Color online.) Left: a single use of the attenuating mem- 
ory channel is described as an elementary transformation which is 
the concatenation of two beam-splitters, respectively characterized 
by the transmissivities /x and k. The first beam-splitter couples the 
memory mode with the local environment, the second one mixes the 
memory mode with the input signal. Right: n uses of the memory 
channel are described as the n-fold concatenation of the elementary 
transformation iB flBll . The concatenation is obtained by identifying, 
for any j, the outgoing memory mode at the j-th channel use with in- 
going memory mode at the (j + l)-th channel use. The memoryless 
limit is obtained when there is no information flow via the memory 
mode, i.e., for /i = 0. 



with ladder operators {aj, a/}j 



I, representing the chan- 



nel inputs, and the corresponding channel outputs associated 
with the ladder operators {&j, 6/}j=i,...n. We model the 
channel environment as a collection of environmental modes 
with {cj, e/}j=i ...„. If the use rate of the channel becomes 
high enough to induce unwanted overlaps between successive 
uses, or to interfere with the finite relaxation time of the chan- 
nel environment, then on a given use the local environment 
could results still "contaminated" by the previous uses. To ac- 
count for such effects, in the same spirit of the Ref.s lul USll . 
we introduce a memory mode that connects all the local envi- 
ronmental modes via an exchange interaction mechanism. 

A single use of the memory channel is modeled by a net- 
work of linear elements: beam-splitters and/or linear ampli- 
fiers being the basic building blocks. Such an elementary 
transformation is depicted in Fig. |2] for the lossy memory 
channel, and in Fig. [3] for the amplifier one lulll . In the 
Schroedinger picture it is formally described by a map 



* {pa,,n) = Tre [[/ {pa^,n ® Pe) U^ 



(3) 



where pa.m is the ingoing state of the input and memory 




FIG. 3: (Color online.) Left: a single use of the amplifying memory 
channel is described as an elementary transfoiTnation which is the 
concatenation of a linear amplifier with gain k > 1 and a beam- 
splitter with transmissivity /j.. Right: the n-fold concatenation of the 
amplifying memory channel is obtained by identifying, for any j, the 
outgoing memory mode at the j-th channel use with ingoing memory 
mode at the (j + l)-th channel use I2l [l3n . 



mode, and p^ denotes the initial state of the local environment. 
The unitary U is the composition of two beam-splitter uni- 
taries for the attenuating channel (Fig.|2]i, or the composition 
of one beam-splitter and one linear amplifier for the amplify- 
ing memory channel (Fig. [3). Finally, the partial trace over the 
local environment accounts for the fact that its final state is ig- 
nored. In the Heisenberg picture, this elementary transforma- 
tion transforms the ladder operators according to 6 = U^aU, 
m' = WmU . For the case of the attenuating memory chan- 
nel, the elementary transformation is the composition of two 
beam-splitters with transmissivities /i, n (transmissivities are 
within the interval [0, 1]), hence we have 



m' = ^J jiK rn + \/l — Ka+ \/(l — ii)Ke , (4a) 

b — s/k, a - a/(1 — ij){l — k) e — ^//(l - k) m , (4b) 

I 



together with the hermitian conjugate relations. Analogously, 
the elementary transformation in the attenuating memory 
channel is the composition of a beam-splitters with transmis- 
sivities /i, and a linear amplifier with gain k, with k e (1, oo). 
The ladder operators hence transform according to 



(5a) 



b = y/na + ^/(l -/i)(K- l)e^ + \/ p.{K - 1) w) , (5b) 



together with the hermitian conjugates. 

The action of the memory channel upon n uses is obtained 
by identifying the outgoing memory mode at each channel use 
with the ingoing memory mode at the following one (see Fig.s 
|2]|3]l. Formally, this is accomplished by assigning amap of the 
form 



^n[pt}j=Tr,[Un 



ra.m 



M 



"»} 



(6) 



J») 



where pa.m denotes the ingoing state of the n channel input 

(n) 

modes and of the memory mode at first channel use, and pe 
is the state of the ri local environments. The unitary oper- 
ator appearing in (|6]l is Un — C/„CAi-i ■ ■ ■ Ui, with Uj be- 
ing the unitary describing the elementary transformation in 
Eq. (|3]l applied at the j-th channel use. In the Heisenberg 
picture the outgoing operators m^^, b/s (for j = 1, . . .n), 
together with their hermitian conjugates, verify the identities 
771^ = Ul^miUn, bj = WajUj. For the attenuating memory 
channel {k ^ 1), by iterating Eq.s (|4|, we have 



i-1 



h=l 



/JIk"' ■' aj + \/{l - p)Ky^^ \fp^ ^ e-j , 



bj ~ y/KUj + y/p, (k — 1) 22 \/a^' 



-j-h-i 



(7) 



j=i 



a-h 



^/{l - p){l - k) ^ ^/pK^ Ch - \/p{l - k)^/pk^ ^ mi . (8) 



h=l 



Similarly, for the linear amplifier memory channel (k > 1), from Eq.s (|5]i we get 



TO, 



^ = -^/Ik" '^1 + Vk - 1 2J VJ^'^ ^ "j + ^(1 - p)k 2J 

■j-h-i 



I n— 7 



(9) 



fej = \/k flj + v^ ('* ^ 1) 2^ vT"^"' 



h=\ 



ah + 1/(1 - m)('^- 1) X! y^^ '^h + ^^(k- 1) V7^^ ^ "ml ■ (10) 
I 



/i=i 



The Heisenberg-picture relations (l7li-(fT0li will be the start- memory channel which will be the aim of the following see- 
ing point of the information- theoretical characterization of the tions . 



Differently from other models where the causal structure 
is not manifest ^ [oL this construction leads to a non- 
anticipatory channel 13] where a given input can only in- 
fluence subsequent channel outputs (i.e. for each j, bj de- 
pends only upon the a^s with h ^ j). The transmissiv- 
ity /i clearly plays the role of a memory parameter. It can 
be related to the ratio between the time delay At between 
to successive channel uses and the typical relaxation time r 
of the channel environment lUOll : for instance we may iden- 
tify /i ~ exp(— At/r). In particular, the model reduces to 
a memoryless (attenuating or amplifying) channel ll22ll for 
/x = (the input aj only influences the output hj), and to a 
channel with perfect memory 111311 for n — 1 (all a^'s inter- 
acts only with the memory mode). These two limiting set- 
tings respectively correspond to the regime At ^ r, and 
At <C T. Intermediate configurations are associated with val- 
ues ^ G (0, 1) and correspond to inter-symbol interference 
channels, for which the previous input states affect the action 
of the channel on the current input [12]. Of particular interest 
is also the case k = where $„ describes a quantum shift 
channel lll2ll . where each input state is replaced by the previ- 
ous one. 

Finally, to exhaustively define the channel model, we have 
to fix the initial state of the local environmental modes 
{ej,eA. Different choices for the environment states lead to 
channels with different features. For instance, an environment 
in a correlated state leads to an additional source of correlated 
noise, similarl y to the effect described by the models consid- 
ered in ]|5l|23l|2^. In the following we assume the local en- 
vironments to be in the vacuum state. With this choice, the 
noise caused by the interaction with the local environmental 
modes is limited to the shot-noise. Moreover, as we will show, 
this choice for the environmental states allows us to unravel 
the correlation in the memory channel jlSll , and to compute 
exactly the quantum and classical cap acity by using known 
results for the memoryless setting 1125112611 . 



Here we only consider the latter case, and procrastinate the 
discussion of this issue to Sec. I VII 

As a first step, we notice that Eq.s (|8]l, ( fTOJ i can be written 
in the following compact form 

bj ^^Ajiiah-^EjiiSh, (k^I), (11a) 

h h 

bj ^J^^jhah + J^Ejhei, (k>1), (lib) 

h h 

where we have defined eg :~ mi, and we have introduced the 
matrices A, E, whose elements are readily obtained from (O, 
([Tol l. In particular, the causal structure of the memory channel 
implies Ajh = Ejh = for j < h. By increasing values 
of n, two sequences of matrices of increasing dimension are 
defined. For each n, we consider the singular value decompo- 
sition of the matrix A, that is. 



ijh 




j'h I 



(12) 



where O, O' are unitary matrices of size n. Actually they can 
be assumed to be real orthogonal, for, in our construction, the 
matrix A has real entries. It follows from Eq.s ( fTTT ) that the 
singular value decomposition of the matrix E reads 



E 



']h 



r=i 



o 



jj' 



(n) 



no]m, 



(13) 



and that 77J"' ^ 1 for k ^ 1, and r/Y^' ^ 1 for k ^ 1. 

We hence define the following set of coUective output vari- 
ables 



bj — Yl '^3'3 ^i' 



(14) 



and the input and environmental collective variables 



III. UNRAVELING THE MEMORY 

The goal of this article is to estimate the classical and quan- 
tum capacity of the model of bosonic Gaussian memory chan- 
nel presented above. As a first step in this direction, we con- 
sider a block of n successive channel uses and show that the 
application of suitably defined encoding and decoding unitary 
transformations allows us to unravel the memory effects, by 
mapping n uses of the memory channel into the direct product 
of n uncorrected channels acting on a suitable set of collec- 
tive variables lUSll . Such a map defines a decomposition into 
'normal modes' of the bosonic memory channel which, under 
certain conditions, is a common feature of bosonic systems 
described by quadratic Hamiltonians (see, e.g., 12711 ). 

Moreover, we ought to distinguish among four different 
notions of (classical and quantum) capacities of the mem- 
ory channel, depending whom the memory mode is assigned 
to ^. Specifically, the initial and final state of the memory 
mode can be under the control of the sender, of the receiver, 
or can be ignored by both and assigned to the environment. 



aj := 5Z^jy" 



i' J 



j:o]r 



(15a) 



(15b) 



These variables are named 'collective' since they are delocal- 
ized over different channel uses. By construction they satisfy 



the canonical commutation relations [b^ , hj 



/] 



Moreover, it follows from Eq.s ( fTTT ) that they 



verify the identities 




(16a) 
(16b) 



We denote Wa, Vb, Te the canonical unitaries 112811 that 
implement the transformations aj -^ a^ — WA aj Wa, bj — >■ 
bj = Vy bj Vb and Cj -^ Cj ~ Te ejTe- We have hence 
shown that the channel <!>„ is unitarily equivalent to the map 



<[pil]=Tr,{z^; 



pin) 
I a.m 



J")' 



u: 



n J ' 



(17) 



a, — 



JfL. 




Wa 


-^a,- 




-^b,- 


Vb 


*b, 

*b„ 


a,^,Kr?''"-^,)^b, 




1 




■^a,- 




-^b,- 


a, J ,^r,«,) *b. 


1 
: ' : 

1 


= i i i 




1 




V7^(^^t/)^b„ 


a«^ -K^^t) *b„ 



FIG. 4: (Color online.) The figure shows a circuital representation of 
Eq. l llSt . Provided that the environmental modes and the initial mem- 
ory mode are in the vacuum state, n uses of the memory channel are 
unitary equivalent to the direct product of n channels independently 
acting on the collective modes. 



„(")' 



J ^(") ' 



with p\ := Te Pe Te, and where the unitary transforma- 
tion U'^ := VsUniWA ® Te) induces the linear transforma- 
tions in Eq.s ( fTSl ). Formally, the unitary equivalence reads 



From Eq.s ^, ( fTOl i. we obtain the following expression: 



M^p) = %' + («,/ - 1) 



//iK 



-\3-3'\ 



where 



S/ 



; + /i(K — 1) 



min{j,i'}- 
/i=0 



{p,K) 



,\h 



(21) 



(22) 



The asymptotic behavior of the sequence of matrices M'"' 
strongly depends upon the value of the product pn. Such 
parameter quantifies the relation between the memory mode 
at two successive channel uses (see Eq.s ( l4al i. dSat). and al- 
lows us to split the parameter region in distinct sectors; on 
one hand, for p.K < 1 the information carried by the memory 
mode is attenuated and 7\/'^"' results in a convergent sequence 
of bounded operators; on the other hand, for pn ^ 1 the influ- 
ence of the memory mode is amplified and the sequence does 
not converge. In the following we will refer to these regions 
as below threshold (fiK < 1), above threshold (pK > 1), and 
at threshold (fiK = 1). 



<[pil] = Vb $„[M^i p^-l Wa] V^ 



(18) 



that is, we can treat the output states of $„ as output of $^ 

in) 

by first counter-rotating the input pa.m by W a (coding trans- 
formation) and then by rotating the output by Vb (decod- 
ing) is]]. Assuming then p^ to be the vacuum state, we have 



Pe 



.(") 



Pe ' and the map (fTTI i can be written as the direct 
product of independent bosonic channels, i.e.. 



$1 = 



^4>{^ 



(")i 



(19) 



J") 



with (\)\rY, ''] being the single-mode transformations whose 
Heisenberg-picture description is given in Eq.s ( fTSI l. A 
schematic representation of the encoding and decoding proce- 
dure which unravels the memory channel is shown in Fig.H] 

In conclusion, we have shown that the n-fold concatenation 
of the memory channel is unitarily equivalent to: n single- 
mode amplifier channels, each with gain m ~^ \ for k ^ 1; 
n single-mode attenuating channels, each with transmissivity 



(n) 



77" ' ^ 1 for K ^ 1; 



IV. ASYMPTOTICAL PROPERTIES 



The aim of this section is to describe the sequence of param- 
eters {jy:' }j=i,... ,„ in the limit of ti — > oo. The regularity of 
the sequence may allow us to determine converging upper and 
lower bounds on the channel capacities, whose evaluation will 
be the subject of the following section. 



For all n, the parameters m 
real, positive matrix 



(") 



are the n eigenvalues of the 



A. Below Threshold 

For /iK < 1 the sequence of matrices A/^"'' is asymptoti- 
cally equivalent ||29:1 to the (infinite) Toeplitz matrix M(°°\ 
whose elements are 



mIT^ := m[°^1 = %, 



jj' 



3-3 



.(oo) _ ^ 



\3-3'\ 



with 



,(°o) 



lim Kjji 

min{jj''}-i-oo 



m(^-i)^ 

1 — pn 



(23) 



(24) 



Following Ref. 112911 the asymptotic distribution of the eigen- 
values of the matrix M^"'^ can then be expressed in terms of 
the continuous function obtained by Fourier transforming the 
elements of the matrix Al''°^\ i.e.. 



77(^)= ^ M 



(°°)gizi/2 



j = -oo 



\/M- V^e' 



z/2 



/pK e' 



z/2 



(25) 



with z €z [0, 2tt] pOQ. The connection between the parameters 



and the function ( 1251 1 is formalized by the Szego theorem 



H29I1 which states that, for any smooth function F, we have 

lini -L^F[^(")]= / ^FHz)]. (26) 



B. Above Threshold 

When the channel operates above threshold the sequence of 
matrices does not converge. For p,K > Iwe find it convenient 
to rewrite Eq. (|2TJ as the sum of two terms: 



M(") := AA^ 



(20) 



Min) ^ g(n)p(«) ^ ^JVf (") , 



(27) 



where P'"^ is a sequence of rank one projectors, c^"^ is a 
diverging sequence of positive real numbers, and AM^") is a 
sequence of matrices which asymptotically converges towards 
the (infinite) Toeplitz matrix AAI^°°\ with entries 



AM 



(oo) _ (l-/i)(K- 1) 



^K 



/^K 



b-f\ 



(28) 



The explicit expressions of the projectors P*^"' and of the ma- 
trices A7\/("^ are reported in Appendix [A] One can easily 
verify that in the asymptotic limit of n — > oo they commute 
(or, to say it more formally, that their commutator is asymp- 
totically equivalent to the null matrix). We conclude that the 
spectrum of the matrices ( 1211 1 is asymptotically composed of 
one diverging eigenvalue [from the diverging sequence c^"-'] 
and of the asymptotic distribution of the eigenvalues of the 
infinite Toeplitz matrix ( |28] l. Similarly to the below thresh- 
old case, the latter can be calculated by Fourier transforming: 
we hence find that it is described by the function in Eq. dZST l, 
extended to the region jik > 1. 



C. At threshold 



At the threshold value, 
expressed as 



/.tK 



= 1, the matrix Af ("^ can be 



M, 



(n) 



%' + (1 - M) + 






min{j,i'}. (29) 



A part from the trivial case /i ~ 1 where M*^"^ coincides 
with the identity operator (perfect channel), the analysis of 
the asymptotic behavior is rather cumbersome (for instance, it 
is not possible to identify a single diverging eigenvalue). We 
thus resort on numerical diagonalization of the sequence of 
matrices. This shows that the finite part of the spectrum is 
again well fitted by the distribution Eq. ( |25l l. extended to the 
manifold with fin ^ 1. 



V. CAPACITIES 

Equation iT% suggests the possibility of computing the 
classical and quantum capacity of the memory channel by 
applyin g the results for the memoryless multi-mode channels 
II25I 12611 . To do so however, we have first to deal with the fact 
that the single-mode channels forming $Jj are not necessar- 
ily identical. Therefore the map iT% is not memoryless in a 
strict sense. To cope with this problem we will construct two 
collections of memoryless channels which upper and lower 
bound the capacity of the memory channel, and then we will 
use the asymptotic distribution (IZST i to show that, for large n, 
they converge toward the same quantity. 

We proceed as follows. For n uses of the memory channel 
we arrange the single-mode channels {'/>(?7,y )} in such a way 

(n) 

that the corresponding parameters rj- are ordered monotoni- 
cally. We further divide the single-mode channels in P blocks. 
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FIG. 5: (Color online.) The figure shows a circuital representation 
of the law of composition of attenuating or amplifying bosonic chan- 
nels. To fix the ideas, let us consider a sequence of n = 3 uses of 
the memory channel, associated to three parameters rjj . For the at- 
tenuating channel, let us assume 771 ^ r;2 ^ 773 ^ 1. It follows that: 
(a) the rate of information transmission cannot be greater than the 
one of the channel with all the transmissivities equal to 773 ; and (b) it 
cannot be smaller than the one with all the transmissivities equal to 
771. The same bounds hold true for the amplifying channel, assuming 
771 ^ 7?2 > 773 ^ 1- 



each containing £ = n/P channels, and define the infimum 
and supremum for each block, i.e. 



(P) 



V. 



vP 



inf 



inf 



n {p-l)e<j^pi 



(n) 



(n) 

:sup sup m 

n {p-l)t<j^pt 



(30) 
(31) 



It is worth noticing that, for any integer P and indepen- 
dently on n, the two collections of parameters, {rjp }p=i....p 
and {7/'^^-'}j,=i^...p, identify two memoryless multi-mode 
bosonic channels 12511 . Such channels are either attenuating 
channels, for k ^ 1, or amplifying channels for k > 1. We 
also remark that (j){rj)(j){rj') = (t>{r]r]'), that is, the composition 
of two attenuating (amplifying) channels with transmissivi- 
ties (gains) 77, ry', is an attenuating (amplifying) channel with 
transmissivity (gain) 77" = 7777'. Since the capacities do not 
increase under composition of channels, we conclude that, for 
any P, the capacity of the memory channel is bounded by the 
capacities of the memoryless multi-mode channels identified 
by the set of parameters {rjp }p=i,...p and {?7''''^''}p=i,...p. 
(See Fig.|5]for an illustrative example.) 

As customary, to avoid unphysical results when discussing 
the classical capacity of bosonic channels, we will adopt suit- 
able constraints on the input encodings (this can be avoided 
when dealing with the quantum capacity, since typically it 
does not diverge even for unbounded inputs i22ll26ll ). Here 
we compute the constrained classical capacity under the con- 
dition 



1 



J=l 



^iV, 



(32) 



where p^"^ is the ensemble of input states over n channel uses, 
the parameter N expresses the maximum number of mean ex- 
citations per mode in average, and the constraint is intended 
to hold for all values of n. It is important to notice that the en- 
coding transformation used for unraveling the memory chan- 
nel [Eq. dlSab l preserves the input energy (this simply follows 
from the fact that the matrix O' is orthogonal). Thus, the form 
of the energy constraint is preserved when written in terms of 
the collective input variables {a^, aU, i.e.. 



iVTr 



.[. 



(")ata. 



^iV. 



(33) 



A. Quantum capacity 

The quantum capacity of the memoryless attenuating and 
amplifying channels has been derived in 12611 . By relying on 
additivity arguments llsill . this result is readily extensible to 
the case of memoryless multimode channels ll25ll . 

For the sake of simplicity, we focus on the limit of un- 
bounded input energy (i.e., A^ -^ oo), which leads to the func- 
tion lEll 



qir]) ■- max{0, log2 77 - logj \v - M} 



(34) 



We hence obtain, for any P, the expressions -p^„q[f]p ] 
and -p ^ qiv^^''] for the quantum capacity (per input mode) 
of the multi-mode channels, respectively characterized by the 
parameters ijp and ri^-^\ Then, we can construct the follow- 
ing upper and lower bounds for the quantum capacity of the 
memory channel: 



){P) 



Q 



{p) 






(k > 1) 
(k> 1) 



(35a) 



(35b) 




FIG. 6: (Color online.) The contour plot shows the quantum capacity 
(measured in qubits per channel use) of the memory channel as func- 
tion of the amplifying/attenuating factor k and of the memory pa- 
rameter /i. The quantum capacity is computed using Eq. \37i . which 
diverges logarithmically f or /i ^> 1 and for k — >■ 1. The dashed lines 
denote the boundaries among different regions, from left to right: at- 
tenuating channel, amplifying channel below threshold, amplifying 
channel above threshold. 



It is worth remarking that the Szego theorem can be applied 
on smooth functions. Since the function in Eq. ( [34b is singular 
for 7/ = 1, we deduce that the expression in Eq. ( l37b coincides 
with the quantum capacity of the memory channel if the pa- 
rameters 1?/^) }, {rip } do not approach the unit value for all 
the values of P and p. By numerical evaluation of the distribu- 
tion of the parameters, we found evidences that this is the case 
for the setting in which both initial and final memory modes 
are assigned to the environment (see Sec.HIIli. As discussed in 
the Sec. IVII the expression in dJTJ i does not coincide with the 
channel capacity if the final memory mode is assigned to the 
receiver. 

Above threshold, the distribution of eigenvalues is com- 
posed by a continuous part, described by Eq. dZSl l. and by one 
diverging eigenvalue. However, for the form of the function 
(I34I 1. the latter does not contribute to the channel capacity. In 
conclusion, the quantity in Eq. dJTl ) expresses the quantum 
capacity of the memory channel both below and above thresh- 
old, proving that Q keeps no record of such discontinuity. A 
plot of Eq. ( [37] i as a function of the parameters jj, k is reported 
in Fig. |6] (for k > 1, the dashed line represents the threshold 
value Kfi = 1). 



so that 



liP) 



i^Q^Q 



{P) 



(36) 



By varying the value of the integer P, we hence obtain a fam- 
ily of lower and upper bounds on the quantum capacity Q of 
the memory channel. 

In the limit P ^- 00, the collections of parameters {t?'^-*}, 

{rjp } may approach a limiting distribution. Clearly, this is 
the case when the channel operates below threshold, with the 
limiting distribution given in Eq. dZSl l. It follows that the upper 
and lower bounds converge towards a same quantity in the 
limit P -^ OD, and by applying the Szego theorem [Eq. ( [26] l1 
we can write 



^TT 



(37) 



B. Classical capacity 

Similarly to the previous case, for the classical capacity 
we consider a family of bounds constructed by introducing a 
proper set of memoryless channels. For P integer, this yields 
the inequalities 



^H 



^C^ c<^^ 



(38) 



where the bounds are computed from the classical capacity of 
the memoryless multi-mode channels defined by the sequence 
of parameters {77'^^)}, {^p '}. 

Let us first consider the case of the attenuating memory 
channel (i.e., k ^ 1) for which exact results can be derived. 
The constrained classical capacity of the memoryless multi- 
mode channel has been derived in |25y, from which we can 
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wnte 



^'"'=pE5feriv,], 



c 



(P) ^ l_ 
p 



p 



Y^a^iP^v]. 



p=i 



where 



g{x) := (x + 1) log2 (a:: + 1) - x logs ^ ■ 



(39a) 
(39b) 

(40) 



The positive parameters {N_}, {Np} describe the optimal 
distributions of the excitation numbers over the collective in- 
put modes. For any P the optimal distributions can be com- 
puted by Lagrange method ll25ll . In terms of a Lagrange mul- 
tiplier L, the optimal distribution is 



^. = {ji: 



^^/< 



-1 



(41) 



and the value of the multiplier is found accordingly to the con- 
straint (33[ . which reads 






Np^N, 



(42) 



and analogously for the distribution {Np}. 

Taking the limit P — )• oo and applying ( |26] l we notice that 
the two bounds (|39] l converge to the same quantity. Therefore 
we conclude that 



C 



dz 

2^ 



g[rj{z)Niz)], (k ^ 1) . 



(43) 



The function N{z) represents the optimal excitation number 
distribution; by taking the limit of Eq.s (HTl i and (02]), it can 
be computed as 



N{z) = [rj{z 



2^/17(2) _ 2 



(44) 



where the value of the Lagrange multiplier is determined by 
the implicit equation 



■Itt 



dz 

2^ 



N{z) ^ N . 



(45) 



In some limiting cases Eq. (03]) admits a close analytical solu- 
tion. For instance in the memoryless configuration /z = 0, we 
get 77(2) = K, N{z) = N and thus coiTecdy C = ginN) Q. 
Vice-versa for k = 1 (noiseless channel) or /i = 1 (perfect 
memory channel) we have ri{z) ~ 1, N{z) ~ N and thus 
C = g{N) (perfect transfer). Finally for k = (quantum shift 
channel) we get 7]{z) = fi, N{z) ~ N and thus C = g{jiN). 
For generic values of the parameters the resulting expression 
can be numerically evaluated, showing an increase of C for 
increasing memory fi. 

Let us now consider the amplifying channel model, ob- 
tained for K > 1. This is intrinsically more complex than 



the previous one, since in this case no exact results are know 
even in the memoryless case. Consequently Eq. (|38l l will only 
provide a lower bound on the real capacity C of our mem- 
ory channel. For this purpose we construct a lower bound 
for the classical capacity C} ' of the memoryless multi-mode 
channel by restricting the coding strategy to only Gaussian 
inputs ll22n (it is worth noticing that such bound is typically 
considered to be tight). This yields 



c(p: 



>^J29[vir\Np + l) + l 



p=l 



v^-i 



(46) 



The optimal distribution of the excitation numbers is com- 
puted by Lagrange method. In terms of a Lagrange multiplier 
L, it reads 



N. 



, = { 



iP) 



1 - 2"-^/'' 



(-P) 



1, 



(47) 



We have hence obtained a family of lower bounds on the 
memory channel capacity. By taking the limit P ^ 00 we 
get 



C^ 



dz 
2^ 



{g{Tl{z)[N{z) + l] + l}-g[i^{z)-l]), 



(48) 
where rj{z) is the asymptotic function ( |25] ), and the optimal 
asymptotic distribution of excitation numbers is given by 



N{z) = [r^{z) 



I _ 2^-^/''(^^ 



- 1. 



(49) 



We stress that, as in the analysis of the quantum capacity Q, 
the expression given in Eq. (08]) holds also above threshold. 
Indeed due to the form of the function (06]l, the diverging 
eigenvalues give a vanishing contribution to the lower bound 
on the classical capacity and can thus be neglected. Figure [T] 
shows the (lower bound on the) classical capacity of the mem- 
ory channel as function of the channel parameters //, k. 



VI. FORGETFULNESS OF THE MEMORY CHANNEL 

According to Ref. 131, when dealing with memory channels 
one has to distinguish among four different settings, which 
may lead to different values of the channel capacities, depend- 
ing whom the initial and final memory modes are assigned 
to. We denote the four possible settings by the label XY: 
XY = AE (the initial memory mode is assigned to the sender 
and the final one to the environment), AB (the initial and final 
memory modes are assigned respectively to the sender and to 
the receiver), EB (initial memory mode to the environment 
and final one to the receiver), and EE (both initial and final 
memory modes are assigned to the environment). If the initial 
memory mode is assigned to the environment, the channel ca- 
pacities may also depend on the memory initialization. [We 
remind that the calculations presented in the previous section 
have been performed in the configuration EE assuming the 
memory mode to be initialized in the vacuum]. 



zi.0.5 




FIG. 7: (Color online.) The contour plot refers to the classical ca- 
pacity (measured in bits per channel use) of the memory channel as 
function of the amplifying/attenuating factor k and of the memory 
parameter /i, under the constraint l l32t with N = 8. For the atten- 
uation channel (k ^ 1), the classical capacity of the memory chan- 
nel is plotted, using Eq. ( 143 1 . For the amplifying channel (k > 1), 
the plotted quantity is a lower bound on the classical capacity. The 
dashed lines denote the boundaries among different regions, from 
left to right: attenuating channel, amplifying channel below thresh- 
old, amplifying channel above threshold. 



The notion of forgetfulness, and hence the class of forget- 
ful memory channels, has been introduced in \2^ for quantum 
channels acting on finite-dimensional Hilbert spaces. Forget- 
ful channels are characterized by several remarkable prop- 
erties, in particular the capacities of those channels do not 
depend on the setting XY and on the initialization of the 
memory mode. Let us notice that our model is defined on 
an infinite-dimensional Hilbert space, and hence the notion of 
forgetfulness cannot be directly applied. However, we can ar- 
gue on the forgetfulness of the channel if an effective cutoff on 
the bosonic Hilbert space is introduced. To prove the forget- 
fulness it is sufficient to show that, in the limit of ti — > oo, the 
final state of the memory mode is independent, in the sense 
specified in Ref. 101, on the memory initialization. Indeed, 
the presence of the exponential factor ,/JIk^ in (|7]), (|9|l sug- 
gests that the channel is forgetful when jin < 1, that is, when 
operating below threshold. We can prove that this holds true 
restricting to Gaussian states with bounded energy. To fix the 
ideas, let us consider the case of the attenuating channel. The 
transformation ^ reads 



/JIk mi + 



X^a 



Y^e. 



(50) 



where a :— (oi, . . . a„)^, e :— (ei, . . . e„)^, and the form of 
the vectors X, Y can be deduced from (|7]i. A Gaussian state 
of the initial memory mode, the input and the environmental 
modes, is characterized by the first moments (7711), (a), (e), 
and by the covariance matrix 



V„, C^ D^ 

V= \ C Fa 

D V^ 



(51) 



where the off-diagonal terms account for possible correlations 
of the initial memory mode with the input and environment 
modes. After n uses of the channel, the state of the final mem- 
ory mode is Gaussian with first moment 

{m'„) = VTIT^" (mi) + XT(a) + Y^ (e) , (52) 



and covariance matrix 

K; = (a*k)" Vm + Va^" (X^C + Y^D + C^X + D'^Y) 

(53) 



X^KX 



Y^TAY . 



The presence of the exponential factors in ( |52] i, ( |53] ) guaran- 
tees the forgetfulness of the channel as long as the energy is 
bounded, e.g., subjected to a constraint of the form 



{m{mi) + X;' 



-.i{a]aj) 



1 



< N. 



(54) 



The same argument can be provided for the amplifying chan- 
nel, as long as it operates below threshold. It is finally worth 
reminding that a first step towards the extension of the no- 
tion of forgetfulness to the domain of continuous variables 
has been made in f23ll . However, we notice that this extension 
is not straightforward if the full infinite dimensional Hilbert 
space is considered. To illustrate this issue we consider the 
computation of the quantum capacity, with unbounded input 
energy, in the generic XY setting. By inspection of Fig.s|2]|3] 
it is immediate to see that in the AB and EB settings, where 
the final memory mode is assigned to the receiver, at least 
one input mode is perfectly transmitted, leading to an infinite 
value of the quantum capacity for any value of the parame- 
ters fi, K. On the other hand, we have shown in Sec. lV Al that 
the quantum capacity is finite in the EE setting (the same can 
be shown in the EB setting) for all values of the parameters, 
provided/^, k 7^ 1. 



VII. CONCLUSIONS 

We have introduced and characterized a model for mem- 
ory effects in attenuation and amplification quantum channels. 
Our findings show that the presence of memory always in- 
creases the quantum capacity of the communication line and 
may increase the classical one. 

Interestingly enough, the highest rates of classical commu- 
nication are reached without the use of entangled codewords. 
For the attenuating channel, the optimal encoding strategy for 
the memoryless channels which bound $„ makes use of co- 
herent states [2511. Since the latter are preserved by the encod- 
ing transformation Wa our results prove, as a byproduct, the 
optimality of coherent state encoding for the presented mem- 
ory channel. 

We emphasis the use of the memory unraveling technique 
II15I1 . which may allows us to evaluate the channel capacity 
without relying on the channel forgetfulness. This shows that 
the unraveling of the memory could have a much broader im- 
pact than the results we have presented here. In particular it 
may find applications on other contexts, e.g., reaching beyond 
current restricted models that involve statistically independent 
errors and are often inapplicable to real physical systems. For 
instance, the model could be easily adapted to describe situa- 
tion in which temporal (causal) correlations between the chan- 
nel uses are replaced by spatial ones or to deal with physical 
models where the underlying noise is more complicated than 
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attenuation (e.g., a stream of two level atoms injected through 
a superconducting cavity). The model can also be applied to 
describe memory effects in quantum repeaters and quantum 
memories with atomic ensembles, where, via the Holstein- 
Primakoff transformation Ii32i1 . the interaction between matter 
and light can be modeled as a formal beam-splitter (or linear- 
amplifier) Hamiltonian luOll . An imperfect swap operation be- 
tween matter and light can be described by a beam-splitter 
with non-unit transmissivity, which in turn may causes mem- 
ory effects. 



i'j 



and 



l^{flK — l) 



AM 



(«) 



(l-At)(K-l) 



flK—l 



flK 



(1 - m) 



k.(k — 1) 



-(1-m) 



2 k(k-1) 



Appendix A: Explicit expressions 

(n) 

Writing the elements of the projector P , as il^jipj'^ the 
explicit expressions for the quantities on the right hand side 
of Eq. ( |27l ) are as follows 



„(") 



k(k-1)(1-jj)^ 

(m«-i)" 



1 



-^^ [(M«)^ 



1 

'1]- 



2(1-A')(k-1) 
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